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EXISTENCE, BLOW-UP AND EXPONENTIAL DECAY OF SOLUTIONS FOR A SYSTEM OF 
NONLINEAR WAVE EQUATIONS WITH DAMPING AND SOURCE TERMS 

LE THI PHUONG NGOC(i), VO ANH KH0A(2), AND NGUYEN THANH L0NG(3) 


Abstract. This paper deals with a system of wave equations in one-dimensional consisting of nonlinear boundary / interior 
damping and nonlinear boundary / interior sources. In particular, our interest lies in the theoretical understanding of the 
existence, finite time blow-up of solutions and their exponential decay. First, two local existence theorems of weak solutions 
are established by applying the Faedo-Galerkin method and standard arguments of density corresponding to the regularity 
of initial conditions. The uniqueness is also obtained in some specific cases. Second, it is proved that any weak solutions 
possessing negative initial energy has the potential to blow up in finite time. Finally, exponential decay estimates for 
the global solution is achieved through the construction of a suitable Lyapunov functional. In order to corroborate our 
theoretical decay, a numerical example is provided here. 


1. Introduction 

This paper is concerned with the following well-known polynomially damped system of wave equations 


(1.1) fmainl} 


utt - Uxx + Ai \ut\'"^ ^ ut = /i (u, v) -I- Fi (x, t ), 
Vtt - Vxx + A 2 \vt\''^~'^Vt = /2 {u,v) + F 2 {x,t) , 


where for i = 1, 2 A^ > 0 the friction terms and > 2 deciding the order of damped parts, are given constants, and fi 
are given interior sources with Fi the external functions which will be specified later. 

Studying systems of wave equations arises naturally within frameworks of material science and physics. Roughly 
speaking, the system (HD is closely related to the Mindlin-Reissiier plate equations (see, e.g., [S]) extending the Kirchhoff- 
Love plate theory for shear deformations, in which they certainly perform three coupled wave and wave-like equations. In 
mathematical aspects, the systems of wave equations have been gaining considerable attention a long time ago and have 
been extensively studying by many authors, see [D11 m ng and references therein. Based on those previous studies, 
many interesting results regarding the existence, regularity and the asymptotic behavior of solutions are obtained. 

Since 2012, Guo et al. [Hj considered the local and global well-posedness of the more general system 


(1.2) leq:sample} 


Utt - Uxx + gi {ut) = fl {u, v) , 
Vtt - Vxx + 52 (Vt) = /2 (u, v) , 


over a connected bounded open domain in with nonlinear Robin boundary conditions on u and zero boundary 
conditions on v. The damping functions gt for i = 1,2 are expected to be continuous monotone increasing graphs 
vanishing at the origin and restrictive in growing up at inhnity. The interior sources impose themselves on the Nemytskii 
operator for the supercritical exponent. Using semi-group equipments, the well-posedness of dm is mostly investigated 
in the relationship of the source and damping terms to the behavior of solutions. One important result included here 
and later in m is that every weak solution blows up in finite time, provides that the initial energy is negative and the 
sources are more dominant than the damping in the system dm. 

In [3J|g, Cavalcanti et al. studied the existence of global solutions and the relationship between the asymptotic behavior 
of the energy and the degenerate wave system/equation with boundary conditions of memory type. The construction of 
a suitable Lyapunov functional, proves that the energy decays exponentially. The same method is also used in [ig to 
study the asymptotic behavior of the solutions to a coupled system having integral convolutions as memory terms. The 
solution of that system decays uniformly in time with rates depending on the speed of decay of the convolutions kernel. 
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In recent years, various types of wave equations with linear or nonlinear damping and sources have successfully solved 
by using Galerkin approximation (see, for instance, [2l[IIl[T3l[T5]). Based on that method associated to a priori estimates, 
weak convergence, and compactness techniques, and via the construction of a suitable Lyapunov functional, the existence, 
regularity, blow-up, and exponential decay estimates of solutions for one of such wave equations is showed in HEIIIS]. 
Moreover, the hnite time blow-up of any weak solutions in which negative initial energy fulfills is obtained in Ca¬ 
in light of the aforementioned works, we put ourselves into the study of the existence, blow-up, and exponential decay 
estimate for the system (ED. 

Let n = (0,1) and Qt = fl x (0,T) for T > 0, a couple of real unknown functions {u,v) is sought for {x,t) G Qt = 
[0,1] X [0,T]. The problem given by (11.11) along with the nonlinear boundary conditions 


(1.3) {eq:main2}- 


u{0,t) =0,-Ux{l,t) + Ki\u{l,t)f^ ^u(l,t) = |ut (l,t)|'^^ ^ut(l,t), 

Vx (0,t) -f K 2 \v (0,t) = fj .2 \vt (0,t)|®“^Ui (0,t) ,u (l,t) = 0, 


and the initial conditions 


(1.4) -[mainS} 


u {x, 0) = uo (x ), ut {x, 0) = ui {x ), 
V {x, 0) = Uo (x ), vt {x, 0) = vi {x ), 


where for i = 1, 2 the constants here Ki > 0, /i^ > 0,pi > 2,qi > 2 are given and for i 
satisfy conditions specified later. 


0,1 the given functions Uj, fy 


2. Preliminaries 

2.1. Abstract settings. Let us first denote the usual functional spaces used in this paper by the notations C™ (fl) , = 

kP™'P (fl) = kpO'P (n) 1 < p < 00 and m G N. Let (-, •) be either the scalar product in or 

the dual pairing of a continuous linear functional and an element of a functional space. The notation H-H stands for the 
norm in and we denote by |1-||^ the norm in the Banach space X. We call X' the dual space of X and denote by 
LP (0, T; A), 1 < p < 00 for the Banach space of the real functions u : (0, T) —>■ A measurable, such that 


and 


\H 


LP{0,T-,X) 



for 1 < p < 00 , 


l|w|lL~( 0 ,r;A) = ess sup ||u(t)|lx, forp = oo. 

0<t<T 

Let u{t),u'{t) = Ut(t) = u{t) ,u" {t) = Unit) = u{t) yVu^t) = Ux{t),Au{t) = Uxx{t) denote u{x,t), ^{x,t), 
^ {x, t), {x, t), respectively. 

On we use the following norm: 

We define 


Vi = {u G : u (0) = 0 } , Y 2 = {v G : V (1) = 0} 

two closed subspaces of . Moreover, the following standard lemmas read the imbedding into (7° (O) and the 
equivalence between two norms, ||ux|| and ||u|fyi on both Vi and V 2 . 

Lemma 1. The imbedding ^ (7° (O) is compact and the following inequality holds 

lbllco(r2) < V^lkllffi > for ally GH^. 

Lemma 2. On Vi and V 2 , two norm v 1 — >■ ||ua;||, and v 1 —^ ||'i'||//i o-xe equivalent. Furthermore, 


||?i|lco(r2) ^ ^ '^1 o-ndY2- 


For the sake of simplicity, we call (P) the problem that contains (11.11) endowed with the conditions (|1.3I) - (I1.4I) . In 
addition, we denote the damping terms and also related functions by {z) = \zf z where r is a given constant. 

2.2. Weak formulation. The weak formulation of the initial-boundary valued problem (P) can be given in the following 
manner: 

Find a pair of real unknown solutions {u,v) belonging to the following functional space 

W={iu,v)€ L°° (0, T; (Vi n x (V2 x P")) : [ut, vt) € L°° (0, P; Vi x V2), {utt, vtt) G L°° (O, T; x P")} , 
such that (u, v) satisfies the variational equations 


(2.1) 


{ e q: Vcir i at i onal< 


{utt {t),4>) + {Ux + Al {ut (t)) ,<(>) + {ut 

= {u (1, t)) (j> (1) -b (/i (m, v),(j)) + (Pi (t), (j)) , 

(vtt (t) , + (vj; (t) +X 2 (^r -2 {vt (t)) , + P 2 d '92 (^'t (0> t)) ^ (0) 


= P 2 ^'p 2 ( 0 , 0 )^( 0 ) + (/2 {U,v),(t>) + (F2 , 


for all 


G Vi X V 2 , together with the initial conditions 


(2.2) {eq: datum} {u (0) ,Ut {Oj) = {uo,Ui) , {v (0) ,Vt (0)) = {vo,Vi) . 

We now make the following assumptions: 

(Al) (uo, Ml) G (Vi n P2) X Vi and (ho, mi) G (V 2 n P^) x V 2 ; 

(A2) Pi,P 2 G pi (0,T;P2) such that F[,F^ G P^ (0,T;P2); 

(A3) there exists P : — )> R the C^-function such that 

dF dF 

(2.3) ■[eq: 2 . 3 } — {u,v) = fi {u,v) , — {u,v) = f2 {u,v) , 

and there also exists the constants a, (3 > 2 and C > 0 such that 

(2.4) F (u, n) < C ^1 -b |m|“ -b , for all u, i; G R; 

(HI) (tio,Mi) G Vi X P^ and (ho, hi) G V 2 x P^; 

(H 2 )Pi,P 2 GP 2 (Qt). 

Remark 3. There are several examples in which functions fi and /2 satisfy assumption (H3), see e.g. dlls]. In particular, 
the authors in [T] considered 


P (tt, u) = q: |u-b -b 2/3 |?«;| ^ , 

where p > 3, a > 1 and /3 > 0. In m, the authors exploited another analytic one 

(2.5) ■[eq:Fcal} p (m, u) = 71 -b -b 72 fyfy , 

where a,/ 3 , 7 i and 72 are positive constants with 72 < 271 . 

3. The existence and uniqueness of a weak solution 

Here we study the existence results involving the uniqueness, of problem (P) which differs from that analyzed by the 
previous works (as showed in the introduction), which consider simpler boundary conditions. It is due to as far as we 
know there are a few papers which take into account complicated types of boundary conditions, e.g. [nma the two-point 
boundary conditions (however, the results for a simple type of nonlinear wave equation are usually carried out). Our 
results included in this section shall therefore extend those previous results with the same strategy, but coupling with 
necessary modifications. 

Now we claim our existence and uniqueness of a weak solution in the following. 



Theorem 4. Suppose that (A1)-(A3) hold and the initial datum obey the compatibility relation 


(3.1) {eq:compal 


I -WOx (1) + (uo (1)) = {ui (1)) : 
(0) + i^2^'p2 ('i'O (0)) = (^'l (0)) • 

If two cases o/(pi,p 2 , 9 i, 92 ) are valid, which are 


\pi,P2>‘2 Ipi,p 2 e {2} U [ 3 , 00 ) 

(2 <91,92 <4, ( 91,92 >4, 

together with the conditions ri,r 2 > 2, then there exists a local weak solution {u,v) of problem (P) such that 


(3.2) {eq :assuml 


' (u, v) e L°° (0, T,; (Vi n p2) X (V 2 n H^)) , 
{ut,Vt)eL^ (0,T,;Vi X V 2 ), 
iutt,Vtt)€L^ {0,T-,L^xL^), 

{QtJ, 


\u' (1, or”' u' (1, •), \v' (0, or”' v' (0, •) G pi (0, n ), 

for T* > 0 small enough. Furthermore, (/91 = 92 = 2 andpi,p2 > 2 , the solution obtained here is unique. 

Proof of Theorem 4. 

Step 1 . The Faedo-Galerkin approximation. Let be a denumerable base of (Vi fl P^) x (V2nP^). The 

approximate solution of (P) is a sequence {(um, in which it has the form 

m m 

'^m (t) = ^ ^ Cjrij (t) 4 >j , Vra (t) = ^ ( dmj (t) , 

1=1 1=1 

where the time-dependent coefficient functions {cmj,dmj) satisfy the following system 

{Um (t) , (fj) + {VUjn (t) , V(t)j) + Ai (d'r-i {Um (t)) , (fj) + {Um (1, t)) (j)j (1) 

= Pidlpj {Ufn (1, t)) (fj (1) + (/l {Um, Vm) i </*j) + {F\ (t) , (f>j) , 

(vm (t) , 4 >j'^ + {yVm (t) + X2 (|d'r-2 (vm (t)) , 4 >j'j + 112^'92 (l^m (0, t)) (0) 

= ^ 24*^2 (I'm (0,t)) (fj (0) + <^/2 {Um,Vm)Aj^ + ^^2 (i) , ^1) , 

^ {um (0) , Ura (0)) = (uq, Ul) , {vm (0) , Vm (0)) = (uq, Ul) , 

for 1 < j < m. A combination of assumptions of this theorem is the direct argument to gain the existence of solution 
Um for the system ( 13 . 31 ) on an interval [OjTm] C [ 0 ,T]. 

Step 2 . A priori estimates. 

Step 2 . 1 . The first estimate. Multiplying the j-th system of ( 13 . 3 L specifically the first equation by Cmj (t) and the second 
equation by dmj (t), then summing up to m with respect to j, and afterwards integrating with respect to the time variable 
from 0 to t, we, after some computations, obtain the following 

Smit) = Smi 0 )+Kif d'pi (Um (1 ,s)) Wm (l,s)ds + P2 / d'pa (I'm (0,s))Um (0,s)ds 

Jo Jo 


(3.3) leq : varil}-< 


+2 


f 


^ Q (Hm (-5) J I'm ('^)) J 11 m (l’)^ ^ (Hm (l') ; I'm (l’)) : I'm (l') 


ds 


(3.4) {eq:Siiill- 
where 


+ 2 f [(Pi (s) , Um (s)) + (P 2 (s) , Vm (s))] ds 
Jo 

= Sm ( 0 ) + Tl + T 2 + 2^3 + I 4 J 








Sm {t) 


iri 

Il^i 

Q 2 


ds 


rt 

= ||wm (0|f + ll^m (t)ll^ + IIVz;^ (i)ll^ + 2Ai / ||um(s) 

Jo 

(3.5) ■[eq:Sm} + 2 A 2 / \\vm ds + 2fj,i / \um {I, ds + 2^12 / 1(0, s)|cis. 

7o io ^0 

By (13.51) and the third equation of (13.31) . there exists a positive constant, says So, that exactly values Sm (0) for all 
m G N, i.e. 

Sm (0) = llUlf + Plf + llVuof + llVuof = 5o. 

We are therefore going to estimate the integrals Ik for fc = 1,4 in the right-hand side of (13.4|) . To do this, let us 
consider the following elementary inequalities. 

Remark 5 (Young-type inequality). Let (5 > 0 and a, & > 0 be arbitrarily real numbers and given q,q' > 1 real constants 
which are Holder conjugates of each other. The following inequality holds 


(3.6) {eq: young} 


ab< + Xs-o'b"^'. 


Remark 6. Given N = ^ max |2; a; /3; then for all s > 0, the inequality s'*" < 1 -I- holds for all 

7G(0,iV]. 

The computations for these kinds of terms li are straightforward. Indeed, estimating Ii and I 2 is the same by using 
two inequalities of the above remarks. Then, by choosing 5 > 0 in such a way that 


5 = 

\ V 2iLi V 2 K 2 

it arises naturally for Ct > 0 representing a bound only depending on T, that 


(3.7) {eq:esll} 


h 


+ Y2 < —Sm {t) + Ct /" [1 + ^m (■^)] ds. 
^ Jo 


Additionally, for Cq > 0 depending only on the initial data uo,vo,ui,vi, and ct,(3, one can prove that 

\Wm it)\\ + \\Vm\\kfi ^ Cq + (a +/3) f [1 -I-5,^ (s)] ds, 

Jo 

the integral I 3 is essentially estimated by the fact that 


I 3 < 2 sup \Tiy,z)\+2Ci+2Ci 

\y\,\^\<Vc^ 

C 

(3.8) ■[eq:esI3} < 


Co + {a + P) [1 -I- Sm (s)] ds 


Co + Co f [1 + S^ (s)] ds, 
Jo 


using also the assumption (A3). The last integral can be found easily by the standard Cauchy-Schwartz inequality: 


(3.9) ■[eq:esI4} 


Y4< lli^ll 


LHQi 


+ 11 ^ 2 ! 


L^Qt) 


f Sm (s) ds < Ct + /" [1 + ^m (®)] 
Jo Jo 


ds. 


We hence obtain from (I3.7I) - (I3.9I) that for 0 < t < 


5. 


{t)<CT + CT[ [l+5^(s)]ds. 
Jo 


We also remark that on the basis of the methods in [9] , there exists a constant T, > 0 depending on T (but independent 
of m) such that 










(3.10) {eq:CT} 


Sm {t) <Ct, Vto G N, VI G [0, T,] . 


This result subsequently allows us to take for all m. 


Step 2.2. The second estimate. We now consider the first equation of (13.31) . Letting t —> O’*", then multiplying the equation 
by Cmj (0) with summing up to m with respect to and importantly by the first compatibility relation (13.1|) vanishing 
the terms [um (1,0)) (1) = (mi (1)) and [um (1,0)) (jij (1) = Ki^p^ (uq (1)) in we thus obtain 

\\um (0)||^ - {Aum (0) ,iim (0)) + Ai |mi Ml, ii™ (0)^ = (/i (uoiflo) , Vm (0)) + (Fi (0) ,Um (0)) . 

Relying on the classical inequalities, we have 


(3.11) ■Ceq:u:} 


Iwm (0)11 < IIAuoll + Ai |mi| 


iri-l 


+ ||/l(MO,ilo)|| + ||i"l(0)||, 


then state that there exists a positive constant, says Ci, that exactly values the right-hand side of (13.111) such that 
\\um (0)|| < Cl for all m G N. 

For the second equation of (13.31) . one also proves without difficulty using similar arguments that there exists C 2 > 0 in 
which it bounds \\vm (0)||) i-e. \\vm (0)|| < C 2 for all m G N. 

Next, we differentiate (13.3p with respect to t, the first equation becomes 


(3.12) 


{'iim {t) , (/)j) -I- (Vlim (t) , -I- Ai (itm (t)) Urn 

= Rl^'pj {Um il,t))Um (l,t) 4 >j (1 ) + (Cr {Um,Vm)u 



and the second one is 


q / -r /X2 {vm it)) Vm (t) + h-2'^'q., {^m (0, t)) Vm (0, t) (0) 

= 1^2^'p 2 (0; ^)) (0) "f i dvdu 


(vm + (VWm {t),V(f>j) + X 

(3.13) ■[eq:3.11} ' ' ' 




for 1 < j < m. Following similar computations above, we multiply the j-th equation of (|3.12|) and (|3.13|) . respectively, 
by Cmj (t) and dmj (t), then summing with respect to j up to m, and in addition integrating with respect to the time 
variable from 0 to t, we shall obtain after some rearrangements 


Fm {t) — Vr 


^ /d'^F , . . , . . d'^F 


(0) + 2^ ( ^ 2 ^m) '^m ('5) 

-a “ ' 


(^mi ^m) ’^rri') ‘^m (■5) / ds 


2 / { — — (Wm; ^m) ('5) “t“ ^ ^m) ('5) ; (5) } ds 


dvdu 


+ 2 / [(^ 1 ' (s) , Um (s)) + (^2 (s) > (s))] ds 


+ 2 / [Rl^'pj {Um (1, s)) Um (1, s) Um (1, s) + K 2 '^'p.^ {Vm (0, s)) Vm (0, s) Vm (0, s)] ds 


(3.14) leq;?!!!!} — Fm (0) F J\F J 2 F JsF 


where 


















Vjn (t) 


11 .. ||2 11 .. ||2 II ||2 II . ||2 

m (0 


+ 


(3.15) {eq;Pm} 


8A]_(n—10 y ^ 

8A2(r2 1) y ^ 

/ ^ (!'“”" I " ^Wm(l,s)^ 

f ^ "" ^^'m(0, s)) 


' 2 

8^1 (gi - 1) '■* 

91 

8Ai2 (g2 -1) 


(is 


ds 


9l 


ds 


ds. 


Combining the arguments from the boundedness of Hiim (0)|| , ||hm (0)|| and the third equation of (13.31) to (13.151) . there 
exists a positive constant, says Vq, that bounds Vm (0) for all m G N, i.e. 


Pm (0) = \\um (0)f + ||h„ (0)f + llVhif + llVDif < Vo- 


(3.16) 

Here the constant Vq is dependent of initial data uq, vq, i(i, Di, the interior sources /i, / 2 , the functions Fi, F 2 and the 
given constants ri, r 2 , Ai, A 2 . 

The technicalities for estimating Vm (t) are almost similar to |13) . so we claim that, by putting 

IC{T,F)= sup \D‘^F{y,z)\, 

|y|,|z|<VCT,|a|=2 

there are possibilities to estimate the first three integrals, i.e. one can show that 


(3.17) {eq:JlJ2} 77i+ 7/2 "F: Ct F j Vm {s) ds^ 

and by the standard Cauchy-Schwartz inequality, we obtain 


(3.18) {eq:J3} J 3 < Ct + [ (|| (s) || + || (s) ||) (s) ds. 

Jo 

For the last integral J' 4 , we shall go through by the following lemma. 

Lemma 7. If one of the following cases is valid, which is 

|pi,P 2>2 |pi,p 2 G {2} U [ 3 , 00 ) 

\2<gi,(72<4, ^ yqi,q2>4, 

the integral J 4 given by can be bounded by 

F 1 

(3.19) {e#tJ4i- 2 / {Um (1, s)) Um (1, s) Mm (1, s) + 7^2 (Vm (0, s)) Vm (0, s) Vm (0, s)] ds < Ct + -Vm (t) . 

Jo z 

Proof of Lemma 7. Since can be divided into two separated integrals for which they can be estimated independently. 
In lieu of taking all possible cases into account, we put myself into the following cases: 

• Case 1: 2 < qi,q 2 < 4,pi = p 2 = 2 and 1 < qi,q 2 < 4,pi,p2 > 2. 

• Case 2: qi,q 2 > I,pi,p 2 > 3 and gi,g 2 > 4,pi = p 2 = 2. 

Now we address below different highly-sophisticated tools corresponding to several definitely complicated terms in each 
case. This endeavor is well accomplished but not rigorously sharp since our strategy is merely to gain a general uniform 
bound (I3.19|) . 

Case 1.1: 2 < qi,q 2 < 4,pi = p 2 = 2 
We have 
























J4 = 2 f [Ki'^2{Umi^,s))Umi^,s)Umi^,s) + K2'^2{Vm{0,s))Vmi0,s)Vmi0,s)]ds 

Jo 

I ll-— I I —-1 

. / \Um (l,s)| ^ Um (1,5) \Um (l,s)| ^ Um (I,-?) ds 

Jo 

1 _ 12. 92 -1 

+ 27^2 / |z^m(0,s)| ^ Vm ( 0 , 5 ) (0,S)| ^ Vm{0,s)ds 

Jo 

f* 2- 

I I UjYi (1, s) I 

Jo 


= 2K, 


AKi 


(3.20) ■[eq:casell} 


9 l Jo 
4 .K 2 


Q ^ I (1 j '5) I (1 ; *5 


92 Jo 




ds. 


Due to (j3.5p and (13.151) which read 


(3.21) {eq:3.20} 


5, 


m{t)>2^if |um (1, s)!*^^ ds + 2/^2 /" | (0, s) |ds, 

Jo Jo 


(3.22) {eq:3.21} 

8/ii (gi - 1) "■* 


Vm {t) > ^Wm(l,s)) 


ds + 


8/i2 (92 - 1) 


9 S Jo 


f 


d / 12. _ 1 

^ (0,s)| ^ Vm (0,S 


ds, 


together with (13.101) . they give us the fact that 


Ji < 


2j^i r 

9i Jo 


^|um(l,s)|'^ '^^+dl ^ (l^tm (1,S)| ^ ^Mm(l,s)^ 


ds 


2i^2 


92 Jo 


/■ 


^|um(0,s)|‘‘ '^''+(52 ^ (I'^m (O^s)! "" %m(0, s)^ 


< 


< 


2 Ki 

qiSi 

2 Ki 

qiSi 


ds 

diKiqi 


[l+|u™(l,s)|^^]ds + ^^ [1 +|u™(0,s)|®^] ds+ 

f SiKiqi 


+ I P™ (i) 


T + - - Sm [t) 

2 fj,i 


2 K 2 

92(52 


T + - - Sm [t) 

2fl2 


- 1) 4y^2 (92 - 1) 

^2^^292 


\4fii{qi-l) 4 /X 2 (92 - 1 ) 


Vm (t) 


(3.23) (.* J4(^ + ^) ft + (4p(‘(0)l' 1) + 4 m!)®i)) (() ■ 


Here we additionally use two classical inequalities which are a'^ '^ < 1 + a'^ for all a > 0, 2 < q < 4 and 2ab < Sa^ + S ^6^ 
for all a, 6 > 0 and d > 0. Hereby, to deduce (|3.19p from (|3.23l) we choose 


dl = ^2 < 


(91 - 1) (92 - 1) 

Kid2qi (92 - 1) + K2^ilq2 (91 - 1)' 


Case 1.2: 2 < gi,g 2 < 4,pi,p2 > 2 
By similar and simple computations, we have 















































Ja = 2 


f {Um (1, S)) ilm (1, S) Ujn (1, s) + K 2 '^'p^ {Vm (0, s)) Vm (0, s) Vm (0, s)] ds 

Jo 

= 2Ki{pi-l)f |Um [wm (l,s)|^ "" Mm, (l,s) (l,s)| "" ^ Um (I, s) ds 

Jo 

n 1 92 -1 

+ 2 X 2 (P 2 - 1 )/ |Mm(0,s)|^^“ |Mm(0,s)| ^ Mm (0, s) | Mm (0, s) | ^ Vm (0, s) ds 

Jo 

= ^Ki{pi-l)J |Mm (1,5)1^"”^ l^m (1,S)|^ "" Mm (1, «) ^ (| Mm (1, s) | ^Mm(l,s))cis 

+ — K2{P2-1)J |Mm(0,s)|^^“ |Mm(0,s)| " Mm (0, s) — (^ | Mm (0, s) | " Mm (0, s) j ds. 

Observe that |Mm (1, s)| + |Mm (0, s)| < VC't by Lemma 2, (j3.5|l and (I3.10|l . one continues by estimating the integral JJa 
that 


J 4 < 


4 £1-1 /"* , ,9-21 d f, ,2l_i \ 

— iLl (pi - 1) C^j? J |Mm(l,s)| " — (^|Mm (1,S)| " Mm(l,s)j 

4 23._1 /"* , ,9-22. 9 i22._i \ 

+ — iL2 (P2 - 1)0,/ J |Mm(0, S)| "" — (^|Mm (0, S)| " Mm (0, s)j 


ds 

ds, 


which gives back to (13.201) in the previous case. By the previously described strategy, (13.191) thus holds. 
Case 2.1: qi,q 2 > 4,pi,p2 > 3 
We start by using integration by parts 


J 4 = 2 (Mm (l,s))Mm (l,s)Mm (1,S) +If24'p2 (Mm (0,s))Mm (0,s)Mm (0,S)] ds 

Jo 

= Ki {pi - 1) |Mm (l,t)|^^“^ |Mm (1,0|^ “ (pi - 2) |mo (1)1^^“^ M? (1) 

pt 

-Ki (pi - 1) (pi -2) |Mm (1, S)|^^“'‘ Um (1, s) Um (1, s) ds 
Jo 

+K 2 {P 2 - 1) |Mm (0,t)|^''“^ |Mm (0,t)|^ “ K 2 {p 2 “ 2) |mo (0)|^''“^Mi (0) 

-1^2 (P2 - 1 ) (P2 - 2 ) [ |Mm (0,s)|^^“'‘Mm (0,s)Mm (0,s)ds. 

Jo 

Still based on (13.51) and (13.101) . we first estimate J4 as follows: 


Ja < K 


< 


{pi-l)SrS -\-Kl{pi-l){pi-2) f 5m" (s) \um j t)f ds 

Jo 

22.-1 I |2 /"* P2-3 , ,Q 

+iL2 (P2 - 1)5„? |Mm(0,t)| + 1^2 (P2 - 1) (P2 - 2) / 5m" (s)|Mm(0,t)| ds 

Jo 

Ot (|Mm (l,t)|^ + |Mm (0,t)|^) + Ot J | Mm (1, s) | ^ + | Mm (0, s) | " 


Secondly, it is immediate to observe that using the inequality < 1 + a'^ for all a > 0 and q >3, together with (13.211) 
and (j3.10|) . let us the following inequality 


rt ^ ^ 

Jo 


ds < T + — ( -1-) Sm (t) < Ct- 

2 \di M2/ 


One also deduces from 






















^Wm(l,0 = |uim(l)h ^ Ulm {^) + J ^ | Mm ( 1 , s) | ^ Mm ( 1 , s)) rfs, 

\V'm{0,t)\'^ Vm{0,t) = \vim {0)\~ Vim (0) -\-J ^(|^m(0,5)|^ Vm{0,s)jds, 

along with the elementary inequality (a + b)^ < 2 (a^ + 6 ^) for all a,b > 0, with Holder’s inequality and (|3.22l) that 


< 2 |Mim (1)1'^" + 2ty (^\um{l,s)\ ^ ^Mm(l,s)^ 


ds 




and by the similar way, we get 


|Mm (0,<)|''' <2|him(0)| 


92 


q^T 


4/i2 (92-1) 


Vm (t) . 


Then, following the natural results of inequalities, which are 

(a + 6 ) ’ < 03 + 69 , Vo, 6 > 0 ,V 9 > 2 , 


we therefore obtain 


ab < 1 — 


S 9-2 a<i-^ 


-S^b^, Va ,6 > 0,V9 > 2,(5 > 0, 

q 


CT(|Mm(l,t)|'+|Mm(0,t)|') < Ct U \u,m (iT + 


qlT 


+Ct ( 2 |Mim (0)1^^ + 

< Ct 


4^1 (91 - 1) 

qlT 


(92 - 1 ) 

2^|Mim(l)|V2^ |hlm(0)r 

o 

2 


Vm {t) 
Vm {t) 


+Ct 

< Co + Ct 

+Ct 


( qlT 


V4^i (gi - 1) 


91 Jl. 

(f) + 


qlT 


4At2 (92 - 1) 


92 _±_ 

Pm" (t) 


9 J 1 UMi(9i-i)y 


91 ( 91 - 2 ) 2 91 ^ 

H-<5]^^ Pm (t) 

9i 


2 \ 32 

1-) 32-2 

92 




V4^2 (92 - 1) 


92 ( 92 - 2 ) 2 92 . 

H- ^ 2 ^ 'Pm {t) 

92 


( JLL \ 

-1- — ) Pm {t) ■ 

9i 92 y 


Hence, to deduce (I3.19|) we choose 5 = (5i = i52 > 0 such that q 26 ^ + qiS ^ 
Case 2.2: qi,q 2 > A,pi = P 2 = 2 

By the same arguments exploited in the previous case, here we can state that 


























J4 = 2 f [Ki'^2{Umi^,s))Umi^,s)Umi^,s) + K2'^2{vTni0,s))Vm{0,s)Vmi0,s)]ds 

Jo 



f^d/\ , ,|2\ 

/■* d , m2 

\ 

Ki 

J jds-\-K2 

f Q CIS V 

j ds 

Ki 

(|m„ (1,7)1"-mL(1)) +K 2 

^I'Mm (0, 7 ) 1 — 

(0)) 

Ki 

1 12 1 1 2 

■Urn (1,7) + 1 X 2 Mm (0,7) , 




also leads to (I3.19|) . 

Hence, we complete the proof of Lemma 7. 

Now, combining (13.171) . (13.181) . and (13.191) we are in a great position to give the fact that 


Vmit) < 2 Vo + 6Ct + 2 [ (l + ||i^{(s)|| + |li^'(s)||)P^(s)ds. 
Jo 

Thanks to GronwalFs inequality, we conclude that 


(3.24) leq: f inalPml T’m (t) < (2Po + dCy) exp 

for all m S N and t G [0, T,]. 


[\l + \\F[{s)\\ + \\Fiis)\\)ds 
Jo 


<Ct, 


Step 3. Passing to the limit. The existence of solution on the interval [0, T*] is now approaching. To summarize, using the 
Banach-Alaoglu theorem (see, e.g., 0), the uniform bounds with respect to m, as stated in the above results p.5l) . (13.101) . 
(I3.15p . and (I3.24|) . imply that one can extract a further subsequence (which we relabel with the index m if necessary) 
such that 


(3.25) {eq:3.24} (um,^^m) ^ (m, v) weak-* in (0,T,; Vi X V 2 ), 

(3.26) {eq;u.v. j(um,'em) ^ (u, u) weakly in (Qt.) x (Qt.) and weak-* in (0,T,;Vi x V 2 ), 

(3.27) -Ceqm..} (itm, hm) ^ (m, weak-* in L°° (O, T*; X L^) , 

(3.28) (um (1, •), (0, •)) ^ (u (1, •), V (0, •)) weakly in (0, T,) x (0, T,), 

(3.29) {eq:3.28} (u^ (1, •) >(0, •)) ^ (u (1> •) > "^ (0, O) weakly in (0, T,) X (0,T,), 

(3.30) (1,-)| "" Urn (1,-): (0r)| " 'I'm (0, •)] ^ (X 1 ,X 2 ) weakly in 77^ (0,T*) X (0,T,), 

(3.31) {eq:3.30}^^ ^ Mm) , ^ (|'i'm|~ ^ ^ (X3,X4) weakly in (Qt.) X (QtJ ■ 

Furthermore, by Aubin-Lions compactness theorem in combination with the imbeddings (0, T») ^ ([0, T»]), 77^ (0, T,) ^ 

(7° ([0, T,]), (0, T^,) ^ (7° ([0, T,]), (0, T,) ^ (7° ([0, T,]), it is straightforward to go on extracting from weak 

convergence results (I3.25|) - (I3.31I) a subsequence {(um,Vm)} such that 

(3.32) ■Ceq:3.31} {um,Vm) {u,v) Strongly in {Qt,) x (Qt,) and almost everywhere in Qt,, 


(3.33) ■[eq:3.32} (^Um,Vm) ^ {u,v) strongly in (Qy^) X (Qt.) and almost everywhere in Qt., 






















(3.34) {eq:3.33} (u^ (1, •),(0, •))(u (1, (0, •)) strongly in (7° ([0, T,]) X C° ([0, T,]), 


(3.35) {eq:3.^^„ (1,-)| " (1, •), [ 1 ;^ (0, •)! ^ (0, •)) (xi> X 2 ) strongly in (7° ([0, T,]) X (7° ([0, T,]). 

Now we have to show the convergence of the nonlinear terms including damping and interior sources. In fact, using 
the continuity argument of /i, one deduces that 

fi —t fi {u,v) almost everywhere in Qt. ■ 

Besides, obverse that ||/i {um,Vm)\\L 2 i^Q^ ) is bounded by ^/^\ sup |/i {y,z)\ which cannot go to infinity, and 

\v\,\z\<VCt 

together with a lemma of Lions [101 Lemma 1.3], one continues to obtain the following 

(3.36) ■[eq;fl} /i (u^, u^) ^/i (m, I') weakly in (Qt.), 

and so is the weak convergence for / 2 , namely 

(3.37) ■[eq;f2} /2 (Mm,Wm) ^ /2 {u,v) weakly in {QtJ ■ 

The weak convergence of damping terms is on the way. Thanks to the inequality 

(zi) - 4-, (Z 2 )| < (r - 1) ki - Z 2 I , Vzi, Z 2 e [-C, C],C>0,r>2 
in accordance with (13.151) . (13.241) and (13.331) . one easily obtains 

(3.38) {eq;psir} (^i^i (um) , (vm)) ^ (^'n (w) , ^'ra (w)) Strongly in (QtJ X {Qt,)- 

It is then worthwhile to mention that (I3.34|) gives 

(3.39) {eq:p4%i {u^ (1, •)), 'Lpa (0, •))) ^ (u (1, •)), ^p, (v (0, •))) Strongly in ([0, T*]) x C° ([0, T*]), 
by the continuity argument of for i = 1,2, and (I3.35|) lets 

(3.40) ■[eq;uv]- (li™ (1, •), (0, •)) Xi, |X 2 |^”^ X 2 ) strongly in (7° ([0, T,]) X (7° ([0, T*]). 

Thus, we take (|3.29l) merging with (I3.40p to have 

(3.41) {eq:3.36} (iXlF”^ Xl, |X2F”^ X 2 ) = (u (1,') , I'(0, •)) 

due to the uniqueness of convergence. 

In the same vein, we claim 

(3.42) {eq:p^®;i^ {'^rn (1,')) : ^ 9 a (0, •))) “t (^gi {u (1,')) , ^92 {v (0, •))) Strongly in C° ([0, T*]) X (7° ([0, T,]) 

from (I3.40|) and (I3.41|) . 

From here on, combining (|3.25l) - (|3.27l) . (I3.32l) - (l3.35p . (I3.36|) . (|3.37[) . (I3.38p . (|3.39l) . (|3.42[) is enough to pass to the limit 
in (|3.3I) to show that {u,v) satisfies the problem (P). In addition, one can use (|3.25I1 - (I3.27|I . (|3.42l) and (A2) to prove 
that 

f Uoox = utt + {ut) - fi {u, v) - Fi £ L°° (O, T,; L^) , 

= Vtt + A24'r-a {vt) - /2 {u,v) - F 2 £ L°° (0,T,;L2) 

which verifies (u, v) £ L°° (O, T ; (Vi D x (V 2 fl H^)) and simultaneously completes the proof of the existence of a 
local weak solution. 












































Step 4- Uniqueness of the solution. The existence of solutions on the interval [0, T,] is established above; therefore, it now 
remains to show the last statement of theorem, namely the uniqueness. Suppose (ui^vi) and (^ 2 ,^ 2 ) are two solutions 
to (P) on the interval [0,T»], which is devoted to the case = <72 = 2 and pi,P 2 > 2, and they go along with the same 
initial datum (ug, ui) and (vg, Di), then they must be equal. We also recall that these solutions must belong to the proved 
functional space properties in p.2l) . 

Define {u, v) := {ui — U2, vi — V2) and based on (12.11) and (12.2L these quantities satisfy the following system of functional 
equations: 


(3.43) ■[eq:3.43} 


{utt {t) , 4>) + {Ux {t) , 4>x) + Ai (d'n (mi (<)) - d'n {u 2 (t)) , (f) + UlUt (1, t) (j) (1) 
= [d'p, (ui (1, t)) - d'pi {U2 (1, t))] cj) (1) + (/l {ui,Vi) - /2 {U2,V2) , (f) , 


Vtt {t) + (vx; {t) ,^x) + A 2 (d'ra {vi {t)) - {v2 (t)) , + pL2Vt (0, t) ^ (0) 


(3.44) ■Ceq;3.44} 

= K 2 [d'p 2 (vi (0, t)) - {V 2 (0, t))] 0 (0) + ( /2 iui,Vi) - /2 iu 2 , V 2 ) , (t> 

for all G Vi X V2 and having zero initial conditions 


w (0) = ^; (0) = Wi (0) =ut (0) =0. 

Taking into account = (ut,vt) in (13.431) and (13.441) . then integrating with respect to t, we obtain the following: 

W (t) = 2 /" {fi{ui,Vi)- fi{u 2 ,V 2 ),Utis))ds + 2 f {f 2 iui,Vi) - f 2 {u 2 ,V 2 ) ,Vtis)) ds 

Jo Jo 

+2Ki f (ui (1, s)) - d/p^ {U 2 (1, s))] Ut (1, s) ds + 2 K 2 [ [dfp^ (wi (0, s)) - d'p^ {v 2 (0, s))] vt (0, s) ds 
Jo Jo 

(3.45) {eq^Wl/Ci + /C 2 + /C 3 + /C 4 , 
where 

W{t) = \\u,it)f + \\vtit)f + \\ux{t)f + \\vAt)f 

+2Ai / (wi (s)) - (u 2 (s)) , ( 5 )) ds + 2//1 / \ut{l,s)f ds 


0 Jo 

(3.46) ■Ceq:Wl} + 2 A 2 [ {^r 2 (^1 ( 5 )) “ {^2 (s)) , Vt (s)) ds + 2/^2 / \vt (0, s)|^ ds. 

Jo Jo 

Our procedure below is essentially similar to the above parts: attempt to estimate ICi for i = 1,4 to derive the 
uniform boundedness of W (t) for which we can use Gronwall’s inequality, then formally make the proof of uniqueness 
self-contained. To handle this, we first state a useful inequality: for all r > 2, there exists Cr > 0 such that 

(d'r- (Zl) - dir (Z2)) (Zl - Z 2 ) > Cr \zi - Z2 T : Vzi, Z2 G K- 
It therefore leads to the fact that 


Wit) > ||«,(t)f+ |k(t)f+ K(t)f+ ||z;.(t)f 

+ 2CrxXi [ \\ut (s)lliVl ds + Cr.,X 2 [ \\vt (5)111% ds 
Jo Jo 


(3.47) ■[eq:WW} 

Secondly, we introduce 


+2fii f |itt(l,s)|^ds-|- 27 t 2 /" |nt(0, s)|" 
Jo Jo 


ds. 


M = max ( II Vw 

i=l,2 




iVt^ill 


L“(0.r.;i?i)) ’ 






z = 1,2. 


Ci (M) = sup 


dh , , 


+ 


dh , . 


Then, it suffices to estimate Ki for i = 1,4. Indeed, applying the Cauchy-Schwartz inequality is straightforward to 
have 


/C1+/C2 <2/ [ll/l (ui,Wi) - /i (U2,U2)|| ||ut (s)|| + 11/2 (ui,Wi) - /2 (U2,U2)|| ||ut (s)||](is 

Jo 

< 2 f (||u (s)|| + ||u (s)||) (Cl (M) \\ut (s)|| + C2 (M) ||ui (s)||) ds 

Jo 

< 2 f (||u, (s)|| + llu,, (s)||) (Cl (M) \\ut (s)|| +C2 (M) ||u* {s)\\)ds 

Jo 

( 3 . 48 ) {eq:KlK2} < 2 (Ci (M) + C2 (M)) f W{s)ds. 


Jo 

To estimate /C 3 and /C 4 , we only need to consider two cases, pi = p 2 = ‘2 and pi,P 2 > 2. First, one may easily show 
that for Pi = p 2 = 2 


/C3 + /C4 = 2iCi 


< 


f u (1, s) Ut {1, s) ds + 2 K 2 f V {0, s) Vt {0, s) ds 
Jo Jo 

2 nt pt J^2 rt pt 

- / il,s)ds + pi / Uf {l,s)ds -\ - - / v^{0,s)ds + p 2 / vf{0,s)ds 

Ml Jo Jo M 2 Jo Jo 


Ki 


(3.49) {eq:K3K4} < [ ^ + ^ ] W (s) ds +-W (t) , 


Ml M 2 / Jo 

where we have followed from (13.471) the intrinsic inequality 


W(t) > (l,t) + (0,t) + 2 ^Mi J \ut{l,s)f ds + p 2 j |ut (0,s)|^ ds^ . 


For pi,p 2 > 2, we have 


/C 3 + /C 4 < 2 


iCi 


(pi-l)C'f^ ^ f |m (l,s)| |ut (l,s)| ds + iC 2 (p 2 - 1) C'r ^ f |u (0,s)| lut (0,s)| ds 

^0 JO 


(3.50) ■Ceq:k3k4} < 


^ (pi - 1) ^ (P2 -1) 


_ Ml M 2 

where we have recalled the following inequality 


1 r* 2 ^ 

I W (s) ds + -W (s), 


(3.51) {eq:3.51} r {.Zi) - ^ r {z 2 )\ < {r - 1) '^\zi-Z 2 \, 'izi,Z 2 &[-C,C],C > Q,r >2. 

Combining (I3.49|) and (j3.50p . we claim that there exists p = p (pi,P 2 ) >0 depending on pi,p 2 such that 


(3.52) {eq;finalK3K4> 


/C 3 + < P (Pl,P 2 ) / W (s) ds + -W {t) . 


Therefore, (I3.48|) and (13.521) together with (13.451) . (|3.46|) and (13.471) imply that 


W(t) <2(2 (Ci(M)+C2(M))+p(pi,P2 )) [ W (s) ds. 

Jo 


Thus, thanks to Gronwall’s inequality, we have W (t) = 0 leading to the uniqueness of solution. Hence, this also 
completes the proof of Theorem 4. 


























Remark 8 . A further fundamental result stemming from the regularity of weak solutions in the above theorem is related 
to the existence of a strong solution. In fact, (lO) allows us to show that there exists a pair of strong solutions {u,v) to 
the problem (P), which satisfies 


' (m, v ) e L°° (0, p*; (Vi n X (V 2 n n c° ([o, r*]; Vi x V 2 ) n ([ 0 , r,]; x l^) , 

{ut,Vt) e (0,r*;Vi X V 2 ) nco ([ 0 ,r*] ;L2 X ^ 

(3.53) {eq: s< G L°° (0,T;L^ X L^) , 


|u'| 2 ^u',\v'\^ ^ v'& {QtJ 


/ {!,■) ,\v' {0,-) G (0,r,), 

Furthermore, if one rules out in some levels the regularity of initial datum, it then also admits a theorem regarding 
the existence and uniqueness of a local weak solution. 


Theorem 9. Taking gi = 52 = 2 and pi,p 2 > 2 into account, let us assume that (A3) still holds. If the initial datum 
which obey the compatibility conditions a only satisfy (HI) and the functions Pi,P 2 o.Te defined in (H2) instead of 
(Al) and (A2), respectively, the problem (P) has a unique local solution (u,v) such that 


r (n, v) e CO ([0, P,]; Vi x V 2 ) n C^ ([0, P,]; P^ x L^) , 

aut,Vt)GH^ (Qt,), 

[uil,-),viO,-)GH^ (0,P,), 

for P, > 0 small enough. 

Proof of Theorem 9. In this proof, let us first recall the weaker assumptions on initial datum and miscellaneous 
functions, that is, (uq, ui) G Vi x L^, (vq, vi) G V 2 x and (Pi, P 2 ) G (Qt) x if (Qt)- Then, we establish sequences 
{(Mom,Mim)} G (fl) , {(z^om, G 0^° (fl) X C(j“ ( ft), and {(Pi^, P 2 m)} G Cg° (Qt) X (Qt) satisfying 


(uom,uim) (uo,ui) strongfy inVi x P^, 
(vom, Vim) (vo, Vi) strongfy in V 2 x P^, 


(Pim, P 2 m) (Pi, P 2 ) strongly in P^ (Qt) x (Qt) ■ 

Note that the sequences {(uom,uim)} and {(wom,I'lm)}, as a result, satisfy themselves the compatibility relation for 
all m G N. So what we can deduce next is the existence of a pair of unique functions (um,Vm) for each m making the 
conditions in the aforementioned theorem self-propelling. Thus, one easily verifies that such functions (um,Vm) for each 
m satisfy the variational problem (I2.ip - (|2.2I1 . namely 


(3.54) {eq :varium} 


' (um (t) , (() -I- (VUm (t) , V(/)) -I- Al (4'ri (um (t)) , 4>) + h-lUm (1, t) (1) 

= PTl^'pi (Um (1, t)) 4) (1) -I- (/i (Um, Vm) , 4>) + (Plm (t) , 4>) , 

(vm (t) , -I- (t) , -I- A 2 (vm (t)) , “k ti 2 Vm (0, t) (j) (0) 

= Pr24'p2 (Vm ( 0 ,t))^( 0 ) + (^f 2 (Um,Vm),^'j + ^P 2 m (t) , , 


for all 


G Vi X V 2 , together with the initial conditions 


(3.55) {eq: inti alum} (um (0) , Vm (0)) — (tlOm, Vim) , (vm (0) , Vm (0)) — (vomt Vim) ■ 

Moreover, the smoothness of (um-,Vm) on the interval [0,P*] is said by (13.531) and we recall below the uniform bound¬ 
edness (independent of m) of Sm (t) on [0, P*] in (13.51) due to the same arguments derived above. 

pt 

Sm(t) = \\um (t)\\^ +\\vm(t)\\^ + \\^Vm(t)\\'^ + \\^Vm(t)\\'^ +2Xi 11 (s) 11 ds 

Jo 

(3.56) {eq:recall} -|-2A2 / \\vm (s)\\^^r 2 ds + 2pi / \um (^, s)\^ ds + 2p2 / li'm (0, s) ds < Cr, 

Jo Jo Jo 







where Ct denotes by a positive constant independent of m and t and t moves along the interval [0,r,]. 
Define Um,k = Um — Uk and Vm,k = Vm — Vk, then these quantities satisfy 


(3.57) {eq< 


for all 


(3.58) {eq:infc5 


Um,k {t) ,4y + {^Um,k (t) + Ai {Um (t)) - (ufe (f)) ,0) + ^ilUm,k (l,t)<^(l) 

= Ki {Um (l,i)) - {uk (l,t))] <))(1) + (/l iUm,Vm) “ fl {uk,Vk), 4 >) + {Fim (t) - Fife {t) , (/)) , 

(ym,k h) + y/Vm,k (t) + X2 {vm (t)) - ^r2 (^'/c (O) i + ^^2Vm,k (0,t)^(0) 

= K 2 [^'p 2 {Vm ( 0 ,t)) - ^p 2 {Vk ( 0 ,t))] 0 ( 0 ) + ^2 {Um,Vm) “ 02 (Wfe,'l'fe) , 0 ^ + (^F 2 m {t) - Fafc (t) , 0 ^ , 
S Vi X V 2 and the initial conditions are 


jU fji 


— i'fkOm ‘^Okj 'aim '^Ifc) , (0) 5 F m (0)^ — ('^Om ^Im ^Ifc) ■ 


It is obvious to obtain the fact that 


Sm,k (t) — Sm,k (0) + 2 f ( '^m) fl i'^k^ '^fc) ; U m,k ^/2 mi ) - 02 {uk,Vk),Vm,k (s)^ 


ds 


+2 


/'[ 


Flm (- 5 ) Fik (s) ^Urn,k (^) / “f \ F2m F2k {s') jFrn,k 


ds 


+ 2Ki / (Um (1, S)) - ^'pi (u/c (1, S))] Um,k (1, s) ds 


+ 2 K 2 / [^'p 2 {Vm (0, S)) - ^'p 2 {vk (0, S))] Vm,k (0, s) ds 
Jo 


(3.59) ■[eq;Smkl} Sm,k (0) + Z\ + Z 2 + Z^ + Z/^, 
where 


^m,k {t') — 


(3.60) ■[eq;Smk} 


Um,k{t) + Vm,k{t) + || VDm.fe (t) || + \\'X7Vm,k ' 
+ 2X1 j yiri {um (s)) - 4 'r'i (w/c (s)) ,Um,k (s)^ ds 
+ 2X2 [ (^r2 (vm (s)) - (vk (s)) , Vm.,k (s)) ds 


+2pi 


Um,k ( 1 ; S) 


ds + 2 p 2 


/■ 


Vm,k (0,s; 


ds, 


(3.61) leq; SmkO}- Sm,k (0) — H'aim 'aifc|| F ||uim “f UVugm ^'aofc|| F UVugm ^'ao/c|| ■ 

The calculations above follow by a valid replacement of the test functions 0 and 0 by Um,k and Vm,k, respectively, 
in (13.571) and then integrating with respect to t. In the similar strategy with classical inequalities such as the Cauchy- 
Schwartz inequality in combination with (I3.51|) . 113.601) and (13.561) . we estimate the terms on the right-hand side of p.59l) 
and obtain 


(3.62) {eq;finalSmk} Sm,k {t) < Cm,k F 2 {1 F r] {pi,p 2 ) F STZt) [ Sm,k{s)ds, 

Jo 

where the appropriately defined terms are given by 


TZt = max sup 


dfi 

dy 


iv^z) 


dfi 

dz 


{y^z) 


V {pi,P2) = —{Pi- 1)^ C't V — {P2 - if ^ 


Pi 


P2 































Cm,k — ‘2Sm,,k (0) + 2 + 2 \\F 2 m — ^ 2 fc|li 2 ('Qy) • 

Here we remark that Cm,k approaches zero as m and k tend to infinity, and by the aid of Gronwall’s inequality, it 
follows from (13.621) that for all t € [0, T»] 

(3.63) {eq;Smklast} (t) < C^^^exp (2r (1 + 77 (pi,p2) + 87 ^t)) ■ 

Thus, one can show that the right-hand side of (I3.63|) shall go to zero as the indexes m and k tend to infinity by the 
direct argument concerning convergences of {(uom,uim)} and {(uom, ^^im)}- Consequently, it gives us the following 

(3.64) ■Ceq:3.64} {um,Vm) ^ {u,v) strongly in (7° ([0,T,]; Vi X V 2 ) H ([0,T,];X L^) , 

(3.65) ->■ (u,v) strongly in (QtJ x (QtJ , 

(3.66) ■[eq:3.66} (u^ (1, •), Um (0, •)) (u (1, ■), V (0, ■)) strongly in (0, T,) x (0,T^). 

One important point should be mentioned here is that by (13.561) we can extract a subsequence of {{um,Vm)} (still 
relabel with the old index m) which reads 

(3.67) ■Ceq:3.67} (um,'em) ^ (m, n) weak-* in (0, T*; Vi X V 2 ), 

(3.68) {eq:3.68} {um,Vm) {u,v) weak-* in (0,T,;L^ x L^) , 

and inheriting from p.64|) - (l3.66p . one deduces 

(3.69) ■Ceq:3.69} {fi {um,Vm), f 2 {um,Vm)) ^ (fi iu,v), f 2 {u,v)) strongly in (QrJ X (Qt.), 

(3.70) ■Ceq:3.70} (^i^i (um) , {vm)) (^'n (w) , 4'r2 {v)) strongly inL^ (Qt.) x (Qt.) ■ 

It is advantageous to state our limit processing now. In fact, pass to the limit in (13.541) associated with (I3.55p . evidenced 
by p.64p - (|3.70p . we obtain a couple of functions {u,v) satisfying the variational problem 

' ^ (“* W > + ('“ 2 = (t) , (l>x) + Al (4'r.i {ut (t j) , (j)) + {ut (1, t)) (j) (1) 

= (u (1. t)) </> (1) + (/i (m, v),cj)) + {Fi (t) , (/)) , 

ft (^) ’ (^) ’ + ^2 (^r2 (^'t {t)) , + M24'92 (Vt (0, t)) ^ (0) 

^ = K2'i>p2 {V (0, t)) ^ (0) -f (^/2 (U, + {F2 {t) , , 

for all ^(/), S Vi X V 2 , endowed with (1^ . of (P) in the current considered case in this theorem such qi — q 2 =2 
and pi,p 2 > 2. This also implies the existence of a local solution claimed by the theorem. The uniqueness of such a weak 
solution is similarly and directly obtained by using the well-known technicalities of regularization procedure investigated 
by Lions (see, e.g., HU)- Hence, we end up with the proof. 

Remark 10. If one has N = ^ max |2; a; /3; | < 1 in Remark 6, and considers the assumptions (HI) and 

(H2), the integral Sm (t) which was uniformly locally bounded by a constant Ct depending only on T (see in (13.101) 1 can 
be estimated globally (and uniformly as well) in time, i.e. 

Sm{t)<CT, VmGN,VtG [0,T],Vr>0. 

Consequently, all arguments used in the proof of Theorem 9 are directly applicable to prove that there exists a global 
weak solution {u, v) to the problem [P) and satisfying 


























(3.71) {eq:regu} 


({u,v)gL°° (0,T;VixV2), 

{ut,vt) € L°° (0,T;L2 X L^) , 

Nevertheless, one point should be noticed is that the case ^ ^ max |2; a; /3; < 1 does not imply 

the weak solution obtained here must belong to (7° ([0,T] ;Vi x V 2 ) fl ([0,T]; x L^) and one cannot also say any 

further statement for the uniqueness in this case. 

4. Finite time blow up 

In principle, blow-up phenomenon of a solution to a time-dependent equation is devoted to the study of maximal time 
domain for which it is defined by a finite length. At the endpoint of that interval, the solution behaves in such a way 
that either it goes to infinity in some specific senses, or it stops being smooth, and so forth. Our main objective here is 
to show that if the initial energy is negative, then every weak solution of (P) blows up in finite time. The result here 
draws from ideas in the treatment of a single wave equation [Tsmi] and also, for example, the recent results for systems 
in mi], but our proofs have to be radically altered. 

Let us first make a brief note concerning the so-called total energy. From mathematical point of view, energy method 
plays a vital role in the study of partial differential equations. That energy is usually used to derive such things as 
existence, uniqueness of the solution, and whether it depends continuously on the data. For a very fundamental scalar 
wave equation which was originally discovered by d’Alembert, one may easily find the energy integral defined by, for 
example, 


i f \utf dx+ i f \Vuf dx, 

^ Jq ^ Jq 

where we basically multiply the equation by ut and integrate over as before using the divergence theorem. 
Physically, the first term in the above energy equation is typically kinetic energy, and the second potential energy. If 
one imposes mixed boundary conditions, says au + bux = 0 on the boundary, where a and b have the same sign, the 
energy can only decrease. If they have opposite sign, the problem is unstable in the sense that the energy will increase. 
In this section, the total energy shall be computed in a common way and then it would be a decreasing function along 
the trajectories, starting from a negative initial value. 

From here on, we aim to consider the problem (P) in a specific case: linear damping = 2 with Fi = 0, and 
Qi = 2,pi > 2, Ki > 0, Ai > 0, Pi > 0 for i = 1,2. Then, we show that the solution of this problem blows up in finite time 
if 

-R (0) = i (lluif + pif + II Vtiof + II VDof) - |fio (l)r^ + ^ |to (^O ix ), UO (^)) dx < 0. 

We also need additional assumptions on interior sources and introduce the total energy associated to the solution (u, v): 
(A3bis) there exists P : —>■ R the C^-function such that 

^ {u,v) = fi{u,v), ^{u,v) = f 2 {u,v), 

and there also exists the constants a, (3 > 2; di, d 2 ,di,d 2 > 0 such that 

(4.1) {eq:4.1} diF {u,v) < ufi {u,v) + vf 2 {u,v) < d 2 F {u,v), for all (u,i;)gR^, 

(4.2) {eq:4.2} di ^|m|“ + |t|^^ < P {u,v) < ^2 Villi'* + |i^|^^ , for all {u,v) G R^. 

(A4) the total quadratic-type energy E (t) is defined by 

(4.3) {eq:energy} 

Eit) = ^ (ll“ wT +11^ wT + wf) - ^ {u{x,t) ,v {x,t)) dx. 




Define H {t) = —E (t) and assume that H (0) > 0, one can show that the time derivative of such a function is never 
negative along the trajectories in a local time, namely 


(4.4) ■[eq:H’} 


H'(t) = Xi \\u {t)\\ +A 2 ||w(t)|| + fJ,i +fi2\v{0,t)\ >0 


for all t G [0,T*) where we have multiplied (11.11) by {u{x^t) ,v{x,t)) and integrated over D. Together with the fact 
that H (0) > 0, we have 


(4.5) {eq:H>0} 0 < H (0) < H (t) , VfG[0,T,). 

Observe the final term in the total energy (|4.3I) . it follows from (14.2|) that 


(4.6) {eq:FF} di (||u(f)||“„ + {u {x,t) ,v {x,t)) < dz (||n(t)||“„ + . 

Combining (14.51) . (14.31) and (14.61) gives 

0<Hit) < (^^|u(l,f)r+^|u(0,f)r)+d2(||n(t)||“„ + |H|^,) 

(4.7) {eq:H<} < (|«(l,t)r^ + \v + ||u(t)||“. + ||u||^,) , 

for all t G [0,T*) and D 2 = max 'fej- It also allows us to derive that for all f G [0,T») 

i(||Vn(t)f+ ||V^(f)f) < (^^|«(l,t)r+ ^|z;(0,t)r)+^J-(«(x,t),u(cr,t)) 

(4.8) {eq:aa} < D 2 (\u + \v + \\u {t)\\lc + \\v\\lfiy 

Now we construct the following functional 

(4.9) {eq:L} L (t) = (t) + e-lp (t) , 

where we define 


(4.10) ■[eq:psiffi (t) := (u (t) , u (t)) + (v (t) , v (t)) 


Y Ik (t)f + Y Ik (t)f + Y“' (k(k^ 


for e > 0 small enough and ^ G ^0, min | | C (O, i). 


To show that one can choose e > 0 small enough such that L (t) is non-decreasing for all t € [0, T,), namely 


(4.11) ■[eq:4.10} L (f) > L (0) > 0, VfG[0,r,), 

a very clear way is to consider the derivative of such a function. In fact, let us prove the following lemma. 
Lemma 11. There exists a positive constant 7 such that 


|P2 


(4.12) it) + ||n (f)f + \\v (t)^ + || Vu (t)r + II Vt; (t)r + Ik WIIl^ + Iklk^ + k (k tT + k ( 0 , t)| 

Proof. First, we multiply HD) by {u {x, t) , V {x, t)) and then integrate over D, the derivative of (t) can be defined as 
follows: 


k' (t) = ||wk)||^ + - (||Vm(<)||^ -b ||Vw(f)f) +Ki |■u(l,<)k' +K2 k(o,t)r' 

+ (/i {u (t) , V (t)) , u (t)) + (/2 (u (t) , V (t)) , V (t)). 


Using this, and the fact from (14.911 that 







L' {t) = (1-0 H-^ {t) H' it) + e (||w (t)f + III; (l)f) - £ (|| Vw it)f + llVi; (t)f) 

+e{Ki +K 2 \v {0,t)f^) + e {{fi {u (t), v (t)), u (t)) + (/z {u{t) ,v (t)) ,v {t))), 

in combination with (14.41) . (14.51) . (14.8p and the following inequality 

(/i {u{t) :vit)) ,u{t)) + (/a {u{t) ,v{t)) > di {u {x,t) ,v {x,t)) dx > didi (||M(t)||^o + , 

shows that 

L'(0 > e (||“ Wf + 11^" Wf) - ^e-Da + |i;(0,t)|^" + ||u(t)||^o + 

£ Ki |w(l,t)|^' +K 2 \v{0,t)f^ +didi (||M(t)||^c + Ikll^p) 

(4.13) {eq:L>} > e (||^, (t) f + ||i; (t) f) + £ (5l - 2 ^ 2 ) (k (1, t)r + (0, t)r + ||li (t)||“„ + ||i;||i,) , 

where we have put Di = min {/4ri,it'2,didi}. If one assumes that 2max 1-^,dal < min {it'i, i^a, we 

deduce 


(4.14) {eq: DSlO < ZI 3 = Di — 2 ZI 2 — min , ^ 2 -, didi | — 2 max 


Pi ’ P2 ’ 



Recall from (|4.7I) that 


min {Ri,R2,didi}. 


D 2 (ill (1, t)r + \v (0, t)r + Ik mi. + mQ >H{t), 

and also thanks to (Ol) and (|4.13p with the assumption to get (|4.14p . it suffices to show that there exists 7 > 0 such 
that (14.121) holds. Hence, we complete the proof. □ 

From this lemma, as said above we obtain (14.111) and realize that the assumption 2 max ■^, da) < min {RTi, RTa, didi} 
is additionally necessary to approach the blow-up result. However, before going to that result, let us consider the following 
supplementary inequalities. 

Lemma 12. Let ^ > 0 such that 

2 < 2 / (1 — 2 ^) < min {a, /3} , 

2 < 2 / (1 - ^) < min{a,/ 3 ,pi,p 2 } , 

then the following inequalities hold 

(4.15) + Ik < 3 (llVll(t)f + |k(i)|li. + k(l,0r^) , 

(4.16) + Ikk)f+ k(0 ,< 3 (iivu {t)f + Ik + \v (0, t)r ), 

for all t G [0, T,). 

Proof. We only need to prove the first inequality since (14.151) and (14.161) are obviously the same. Put si = 2/ (1 — 2^) 
and sa = 2/ (1 — ^), our strategy is to independently consider each terms on the left-hand side of (14.151) and furthermore, 
to estimate those quantities we mainly divide into two cases. For clarity, let us first consider two cases for |k|lic. 

• Iklli" ^ 1 - Oil® easily deduces from 2 < si < a that 

Ik wiIl. < Ik ml. < II vii (t)f + Ik it)\\i. + k ( 1 , t)r . 

• lkllL“ — Similarly, we get 


Vw G Vi, 
Vi; G Va, 


















ih m% < Ik < IIwf + Ik wii^ + k (1, tT • 

Next, we consider two cases for ||u||: 

• ||m|| < 1: Since 2 < S 2 < a, 

ikwir^ < ikk)f < iivu(t)f+ ikk)iii. + k(i,i)r • 

• Ikll — 1- III III® same manner, we have 

Ik {t)\r < Ik wr < Ik wiil. < II vu k)f + Ik wiil. + k (i,i)r • 

Finally, we go through the last term |m (l,t)|: 

• k(l,I)l < 1: 

|u(i,t)|*" < \u{i,t)\^ < lk(i)llco(n) ^ llVuWf < IIVM(t)f + IkWIli- + k(i:i)l^' ■ 

• \u (1, t)| > 1: It is natural to say \u (1, t)k < k (1) 1)1^^- 
Combining all of the above inequalities completes the proof. 

□ 

Theorem 13. Suppose that (A3bis) holds and H {0) > 0. //2 max | J2I < min {itTi, 1^2, <^ik} in (A3bis), 

then any weak solution {u,v) of (P) blows up in finite time for any (uo,ui) and (hojUi) in (HI). 

Proof. By using an elementary inequality 

Vr > l,a:i > 0,i =1/7, 

\z=l / i=l 

and thanks to (|4.9p . ()4.10p . one can show that 

(t) < C (^H (t) + \{u{t) ,u{t))\^^^^ + |(u (t) ,h 

(4.17) ■[eq:Lxi} + ||u + |k (1,1) + (0,t)) , 

where C = max , (p2£)^^^^~'^^|- 

Moreover, we find 

(4.18) ■Ceq:4.17} | (u (t) , U (t)) |< IkWir^^'”^^ || U (t) f < |k W ^ W , 

where we have used the Cauchy-Schwartz inequality. Holder’s inequality. 

According to Young’s inequality introduced in ()3.6L by choosing 5 = \,q = , q' = 2(1 — ^) and letting a = 

Ik (I)IIl“^^~^^ k = ||ii (k||^^*'^ 11^® easily obtains from (|4.18l) that 

(4.19) ■Ceq:4.19} | (u (t) , U (t)) |< Ci ^|k (t)||^'i‘'^“^^^ + ||u (t)||^^ , 

where ci = max | 2(i~-|) ’ 2(1-^) } ^ I- Similarly, we have for a constant C2 £ (0,1) 

(4.20) {eq:4.20} \{v (t) ,V {t))\^^^^ < C2 (|k (l) + ||i’(0 |f) • 

Therefore, the observability estimate (14.171) along with (I4.19F (I4.20|) . (I4.15p . and (14.161) implies that there always exists 
a positive constant (here we reuse the notation C for simplicity) such that 

L^/P-i)(t) < C(p(t) + ||u(t)f+ ||u(t)f+ ||Vu(t)f+ ||Vu(t)f 

+ Ik WII^ + Ik k)ll^. + k (1, t)r + k (0, t)r) , Vt G [0, n) . 


(4.21) {eq;LL} 



















In addition, the estimate (14.1211 together with (14.211) allows us to take a positive constant C such that 


(4.22) {eq:4.22} L'(t) > (t) , Vte[ 0 ,r,). 

Now, integrating (14.221) over (0,t), one has 


Li/i^-0 (t) > 


t e 


L-«/(i-«) (0)- 

which yields that L (t) definitely blows up in a finite time given by 

l-Cr- 




(0) , 


T, = 




-L-«/(i-«) ( 0 ). 


Hence, we complete the proof of the theorem. 


□ 


5. Exponential decay 

While the previous section employs the total energy to find the finite time which produces the blow-up result, the total 
energy in this section is considered to be in exponential decay. In particular, we study the global solution (u, v) of (P) 
satisfying (13.711) and corresponding to ri = r 2 = 91 = <72 = 2 and pi,p 2 > 2, and find an related-initial data in the stable 
set together with suitable and necessary assumptions provided by, for example, external functions, for which the solution 
decays exponentially. Like the blow-up phenomenon, this sort of results can be seen as an extension of many previous 
works from the single wave equation in, for example, mm to the system of equations. 

Our result here relies on the construction of a Lyapunov functional by performing a suitable modification of the energy. 
To this end, for 5 > 0, to be chosen later, we define 


(5.1) ■[eq:Lcal} £ (t) = E {t) + Sip (t) , 

where we have recalled the energy in ( A4) and the function ip (t) in (14.171) . 

Under some additional conditions, it suffices to see that C{t) and E{t) are equivalent in the sense that there exist two 
positive constants /3i and /32 depending on d such that for t > 0, 


(5.2) {eq: equi} PiE (t) < /I (t) < P 2 E {t) . 

Before explicitly providing this equivalence, let us first consider the time derivative of the total energy, which can be 
defined in a similar way for (14.41) by the following lemma. 

Lemma 14. The time derivative of the total energy satisfies 


(5.3) {eq:E’ll 


E'(t) < —A* ^||m (t)||^-I-||z; (t)||^^ — p* (l,t)|^-I-[h (0,t)|^^ 

I m + \\f2 (t)ii) + i (iiPi m + \\e2 (t)ii) (||u (t)f + \\vit)f) , 


E'(t) < - ^||rt (t)||^-I-||r; (t)||^^ - p, (l,t)|^-I-|r; (0,f)|^^ 

(5.4) {eq:E>2> (||Fi (t)f + ||P 2 (t)f) , 

for all El > 0, A* = min {Ai, A 2 } > 0, and p* = min {^ 1 ,^ 2 } > 0. 

Proof. Let us remark that multiplying (ED by {u{x,t) ,v{x,)), integrating over H and using integration by parts, we 
obtain 


(5.5) E' (t) = -Ai ||u(t)||^ - A 2 ||t (t)||^ - Pi |u(l,<)|^ - P 2 |p (0,t)|^ + {El (t) ,u{t)) + {F 2 (t) ,v{t)) . 

To obtain (15.3p . we only need to estimate the last two terms. It is straightforward that by using the standard inequalities 
which read 
















{F2{t),v{t)) < ^\\F,it)\\ + ^\\F2{t)\\\\v{t)\\\ 

it yields 

(Fi + {F, i m (i)ii + m wii) + i m wii + ii^’2 wid (of + Ht; wf). 

For the second estimate (15.41) , we also use the same approach to prove. Indeed, by Young’s inequality one deduces the 
following 


{F, {t),u{t)) + {F, {t),v{t)) < ^ (||Fi {t)f + \\F, (t)f) + I (||u(t)f + ||r; (t)f) . 

Hence, the lemma is proved completely. □ 

Let us next define the following functions L (t) = li {u (t)) for i = 1, 2 and J {t) = J {u (t)) by rewriting the expression 
of the total energy E (t) in (14.31) : 


(5.6) 

■[eq:Enew} 


Eit) = ^[\\u{t)f + \\v{t)f)+J{t), 

(5.7) 

feq:J} 


^ W = J f 1 - - - (11 Vu (t)f + II Vu (t)f) + ^-^ + 

2 \ Pi P2 J \ 7 Pi P2 

(5.8) 

{eq:11} 

h{t) 

= \ (l|VM(t)f + liVu {f)f^ - Ki\u{l,t)f^ - Y J^^i'>J'{x,t) ,v{x,t))dx, 

(5.9) 

■[eq:I2} 

hit) 

= ^ (11 Vm (t)f -f II Vu (t)f) - K 2 \v ( 0 , ^)r" - y ^ ^)) 


Furthermore, let us provide a further assumption as follows: 

(H2bis) Fi,F 2 G (iR+iL^) (K+;L2). 

From now on, our main result in this section is established where the proof of the equivalence between C {t) and 
E (t) is also included. It says that if there is an exponential rate of energy decay for external functions Fi and F 2 , and 
influenced by such functions, the related-initial energy function, says i?,, is properly suited in a particular set, then the 
quadratic-type total energy decays exponentially. 

Theorem 15. Suppose that (Al), (A3bis) and (H2bis) hold, together with the d 2 < min{pi,p 2 } in (ASbis). Assume 
that Ii (0 ), I 2 (0) > 0, and the energy given hy i4-^ such that its initial energy satisfies 


(5.10) {eq:eta*} ry* = ^-^——^d 2 


{p^E^)'^ ^-I--I-IFi 2 ^+K2{p»E^)'^ ^ < 1, 


where p* = , E^ = [E (0) -I- p\ exp{2p) and p = \ (HF’i (s)|| + \\E 2 (s)||)(is. Moreover, let us assume the 

external functions decays exponentially in the sense that 


(5.11) {eq:FlF2} ||Fi (t)||^-f IIF 2 (t)||^ < ?yi exp (- 772 ^) , Vt > 0, 

where pi and p 2 ore two positive constants. Then there exist positive constants C and 7 such that 


(5.12) {eq: final! 


E (t) < Cexp {—"ft ), Vt > 0. 








Proof. Firstly, we claim that p (t) > 0 for i = 1,2 and for all t > 0. In fact, since li {t) ,i = 1,2 is continuous and its 
initial value is positive, thus there exist two positive constants Ti and T 2 such that 

h (t) >0, Vt e [0, Ti], and h (t) >0, Vt e [0, T 2 ] 
which also leads to the fact that 

(||Vu(t)f+ ||Vu(t)f), vte [o,r], 

where T = min{Ti, T 2 } or we can say that 

(5.13) {eq:5.10} ||Vu(t)f + ||Vu(t)f < -- J (t) < -- E (t) , Vter0,rl. 

P 1 P 2 - P 1 -P 2 P 1 P 2 - P 1 -P 2 

By (15.31) and thanks to p (0) > 0 for i = 1,2 associated with (j5.6L we deduce 

E{t)<Ei0) + -J^ {\\F,{s)\\ + \\F 2 {s)\\)ds + m{s)\\ + \\F 2 {s)\\)E{s)ds, Vte[0,T], 

then using Gronwall’s inequality leads to the following 


E{t) < 


E{0) 


1 r 
sio 


l^^lWII + ll^’2(s)||)ds 


exp 


|i^lW||+||i"2(s)||)ds 


(5.14) {eq:5.11} < (0) + p] exp (2p) , 

where we have defined p = 5 /q°° (HF^i (s)|| + 11^2 (s)||) ds. 

Therefore, we combine (15.141) with (|5.13p . and put i?* = [E (0) + p] exp (2p) ,p* = to obtain 


(5.15) {eq:5.13} || Vu (t) f + || Vu (t)f < , Vt G [O, T] . 

Observe the assumption (14.21) in (A3bis), and also (I5.15p . it yields 


K,\u{i,t)r+K2\v{o,t)r 

(5.16) {eq:5.14} 

and 


< Ki\\vuit)r+K2\\vv{t)r 

PI -I rj P2 1 O 

< Ki (p^E,)--^ II Vu (t)f + K2 (p,f;,)--^ II Vu (t)f , 


{pi+P 2 )J T {u{x,f) ,v{x,t))dx < {pi + P 2 ) d 2 [\\u {f)\\- 

(5.17) {eq;5.15} ^ ~ ^ 

Thus, it follows from (15.161) and (15.171) that 


< {pi+P2)d2 (p*A,) 2 “-^ ||Vu (t)ll^ + (p,F;,) 2“^ ||Vu (t)ll^ 


Ki |m (l,t)|^^ + K 2 \v (0,t)|^^ < P* (||Vu(t)||^ + ||Vu(t)||^^ 

< (||vu(t)f+ ||vu(t)f), vtG[o,r], 

where p* is given by (15.101) . 

So we now claim that Ii (t) and I 2 (t) are positive for all t G [O, T]. If we put 


T, = sup {T > 0 : /i (t) and I 2 (t) are positive Vt G [0, T]} , 

and if T* < 00 , then (by the continuity of Ii (t) and I 2 (t)) we have Ii (T*) and I 2 (T*) are non-negative. Therefore, by 
the same arguments, it is possible to show that there exists f* > T* such that p (t) > 0 for i = 1,2 for all t G [O,!).]. 
Hence, we can conclude that p (t) >0 (i = 1, 2) for all t > 0. 

We shall next end up with the equivalence of £ (t) and E (t). Let us recall that by (15.11) . (j5.6l) . (|5.3I) . (15.4p . (j5.7p and 
(14.101) the expression of £ (t) can be rewritten as follows: 
























h (t) , h {t) 


Cit) = I (||« (t)f + lit; Wf) + i (^1 - 1 - {\\Vu {t)f + II Vz; (i)f) + ^ + 


JAi II /,\\\2 *^^2 


+S{u(t) ,u(t))+S{vit) ,vit)) + — ||u(t)|| 


P2 


It is straightforward to see that 

1 


C{t) < i(l + 5)(||u(t)f+ + 

(5.18) {eq;upper} +— - - -h (5 (Ai + A 2 + + ^ 2 )^ (^)ll^ + 11^^ (^)ll^) ) 


P2 


where we have used an elementary inequality 

1 


{u{t),u{t)) + {v{t),v{t)) < - (||Vu(t)f + ||Vu(t)f) + - (11 Vu {t)f + II Vu (t)f ) . 
Therefore, one can choose from (15.181) that 

1 ~ ^ ^ ^ (-^1 + -^2 + Ml + M2 ) 


(32 = max 11 + 5, 

to obtain C (t) < (32 E (t). 
Furthermore, one easily finds that 


pi P2 


= max < 1 + (5 ,1 + 


i5 (1 + Ai + A 2 + Ml T M2) 


1 — - — 

pi P2 


cm>l(i-S )(||„(,)||%||.(*)||^) + A + A + 1 (1 _ 1 _ 1 _ (11 v„(,)f + Iiv.(Of), 


Thus, choosing <5 > 0 small enough for which 

(3i = min < 1 — 5 ,1 — 


1 _ A_ A 

Pi P2 


> 0 , 


most likely we choose S G I 0,1 — j°_ j_ J to obtain £ (t) < I 32 E (t). Hereby, our Lyapunov functional £ (t) is 

\ PI P2 / 

definitely equivalent to the total energy E (t) for alH > 0. 

It remains to consider the functional ip (t) which is defined by (j4.10|l . The time derivative of ip (t) can be found by 
multiplying CH) by {u (a:, t ), v (x, t)) and then integrating over H. It therefore has the following expression 

Ip' it) = ||zt(t)||^ + ||^'(^)||^ - l|VM(t)f - ||Vz;(t)||^ + iLi |M(l,t)|^" +1^2 |u(0,t)|*’'' 

(5.19) ■Ceq;5.17} + (/i {u (t) , v (t)) , U (t)) + (/2 (u (t) , v (t)) , V (t)) + (Fi {t),u{t)) + (F 2 (t) ,V (t)) . 

On the one hand, we have 

if I {u (t) ,v{t)),u (t)) + (/2 {u (t), V {t )), V (t)) <d 2 [ E{u (x, t ), V (x, t)) dx, 

Jn 

and by (15.81) and (15.91) . we continue to estimate the above inequality 


{fi(u{t),v(t)),u{t)) + {f 2 {u{t),v{t)),v(t)) < d 2 


U- + -) (llVix(t)f+ ||Vz;(t)|r 

2 V Pi P 2 / ^ 


(5.20) {eq:5.18} 

On the other hand, one easily has 


-(^\uii,t)r + ^\vio,t)r]- 

, Pi P2 


h (t) , h (t) 


Pi 


P2 


(5.21) {eq: ^ili9(}t) ,u(t)) + (£2 (t) , v (<)) < y 


Vu(t)f + llVz; (t)f) + ^ (||Fi (t)f + IIF 2 (t)f) , V£2 > 0. 


















Thus, we obtain from (I5.19II - (I5.21|) that 


it) < Ih wf + llWOir - (l - I - y (^ + ^)) (l|Vu(t)f + ||Vz;(t)f) + (l - ^) Ki 

(5.22) ■Ceq:estpbi^}l - \v (0,1)^ - ^2 ^ (H^ {t)f + WF^ (t)f) . 

Now the complete proof is coming. By (15.41) . (15.221) . and thanks to (15.161) and (15.101) . the time derivative of our 
Lyapunov functional C (t) can be estimated as follows: 


£'(t) < - 


,t)\ + |u(0,t)| j + — 


(a* - y) + ||f (t)||^) - /i* (^|M(l,i 

+A (||u(t)f + “^(^“y“y(y + (ll^^ +11^^'Wf) 

/ j \ 1 


^K,\u{l,t)r + (^l-^'^K2\v{0,t)\ 


< — ^A* — S — 


■) (ih(^)f + 11 ^Wf) “ 


+5 

-Sd2 


(^1 - K, ||Vu (t)f + (^1 - K2 IIVu (t)ll' ' 

/ r /■_ /'+^\ 1 / 1 ;c \ / _ _v 


h jt) _|_ h jt) 

Pi P2 


1 / j_ 

2 \ei £2 


< - 


(||y(t)f+ ||£2(t)|| 

\ t'i- '-'^Z ^ 

(A*-5-y) (||u(t)f+ ||u(t)f) -^2 




(||y(t)f+ ||F2(t)ir 


(5.23) {eq;L’eaS} (1 — 77») ( 1--- 

V max{pi,p2| 

for all i5, £1,62 > 0. Here we imply 


1] (iivu(t)f+ iivu(t)ir 


^2 < min{pi,p2} , 0 < £2 < 2(1 - ?7,) ( 1- p- -^ . 

V max{pi,p2|/ 

Then for 5 small enough, satisfying 0 < (5 < A*, and let £i > 0 such that 0 < £i < 2 (A, — 5), the equivalence of C {t) 
and E {t) together with (|5.23l) and (15.111) says that there exists a constant 7 £ (0, 72) such that 

C' (t) < -7/3 (t) + fyexp {-p 2 t ), Vt > 0, 
and by Gronwall’s inequality, we have 

£ {t) < ^£ (0) + —^ exp (-yt), Vt > 0, 

which leads to (|5.12|) by the equivalence of £ {t) and E (t). Hence, we complete the proof of the theorem. □ 

A numerical example. The one-dimensional linear damped system of nonlinear wave equations has been qualitatively 
investigated in the sense of exponential decays and therefore, in this subsection the emphasis is put on the illustrative 
framework. We show below a numerical example where we recall that our considered problem in this section is particularly 
given by 


Utt - Ua;x + Alltt = fi (it, v) + El {x, t) , 
Vtt - Vxx + A2Ut = /2 (it, v) + F 2 (x, t) , 


































along with the boundary conditions m and initial conditions (m. We notice that the interior sources /i and /2 
have also been indirectly introduced in Remark 3. In fact, the functional G (K^; K) which is observed in (12.51) as an 
example, can be established by 


fi{u,v) = a (71 l^r ^ + y |u| 2 
f2{u,v) = 

By these explicit expressions, (lOl) obviously holds. Moreover, we deduce the facts that for all (m, v ) G 

(71 - y) (kr + (u,v) < (71 + y) (l^r + , 

min {a, /?} T {u, v) < ufi {u, v) + vf 2 (u, v) < max {a, /?} T {u, v ), 
which imply di = min {a, /3} , ^2 = max {a, /?} , di = 7 i — and ^2 = 7 i + ^ in (14.11) . (|4.2I) and if one can choose an 
appropriate set of constants a,/ 3 , 7 i, and 72 such that 72 < 271 and d 2 < min{pi,p 2}5 (A3bis) is clearly valid. 

In our example, we take a = /3 = 4, pi = p 2 = 6 and Ki = = Xi = 1 for z = 1, 2 with 71 = |, 72 = 5 (so, d 2 = i and 

(^2 = !)■ Next, the initial conditions are defined by 

uq (x) = x (e® + 1 ) , ui {x) = —xe® (e® + l) , 

(a;-l)e® (e® + l)^. 

{x — 1) (4x^ — 6a; + 3 ) 56®®“** (1 — a;) 

(e9+4z + 1)3/4 ^^ 9/4 ■ 

(5.24) {eq:exact} We. (x, t) = ^=^_, v^^X.t) = 

but they shall be neglected since we want to consider the illustrative approximation solutions. 

In order to fulfill all assumptions of Theorem 15, it remains to check (15.101) and (15.111) . Obviously, one has p* = 3 and 
p is computed by the following 


ho (a;) = (1 - a;) (e® + l) 
Then, our external functions are 


hi (x) = 


Fi (x, t) = - 


4x3 — 2 x® + : 


5e®+^*a 


.9+44 , 


1 ) 


3/4 


(e9+4t + 1) 


9/45 -^2 (a;, t) 


Finally, the exact solutions here are given by 


P = 


< 


1/°°(i+.(»)ii + iift(»)ii)* 

f 
f 

Jo 


25e2(9+4.) 


19g®+4^ 


105 (e®+4* + 1)^/^ 3 (e®+4* + 1)®/® 3 (e9+4* + 1)^ 


1/2 


ds 


-|(9+4s) , ^ |(9+4s) , l^g-2(9+4s) 

105 3 3 


1/2 


ds 


/■°° /168 2 27 g\ 

(5.25) 4.„6/25) < 


1/2 


ds < 1.563 X 10 


-3 


and the initial energy can be estimated as follows: 


E{0) 


^ +ilof+iivhof+iivhof) - (iho (i)r+1^0 (o)r) 

~i/ ( 1^0 ( 2 ^) 1 ^ + 1^0 (a;)|'‘) dx-iy |ho(x)|^|ho(x)|^dx 
lgi 8 (g 9 + i)-i + (g9 + 1 )^ _ 2 (e9 + (e® + 1 )-' - 

o J-U 


1 


(e® + l) ^<0.015. 


(5.26) {eq;l/4} 




















Combining (15.251) and (15.261) . the indirect argument of checking (15.111) is made by (15.251) . and we obtain that 


i5» = [E (0) + p] exp (2p) < 0.017, 

which leads to 


ry* = 12p»i?* + 2 < 1. 

Therefore, the assumptions needed to check are all satisfied. 

At the discretization level for this problem, a uniform grid of mesh-points (a;^, t„) is used. Here Xk = kAx and tn = nAt 
where k and n are integers and Aa; = , At = ^ the equivalent mesh-widths in space x and time t, respectively. We 

shall first consider the following differential system for the unknowns {Uk (t), 14 (t)) = {u {xk,t ), v {xk, t)) for k = 0,K: 

(5.27) feq: 5. 26} ^ (t) = {U 2 (t) - 2t7i (t)) - Hi (t) -f 3t7f (t) -b (t) Ui (t) -b El {xi , t), 

(5.28) ^ (t) = (t/fc_i (t) - 2t/fc (t) + Hfe+i (t)) - Wfc (t) + 3t7| (t) + 14 ^ (t) (t) + (a:^, t), fc = 2,/f-l, 

(5.29) ^(t) = K {Ul, (t) - Wk (t)) - if" (Hk (t) - (t)) - Uk (t) + 3[/i (t) + i^i (a^^^), 

and 

(5.30) ^ it) = K {V^ (t) - Vo (t)) + if" (Vi it) - Vo (t)) - Vo (t) + 3Vo3 (t) + F 2 (a;o, t), 

(5.31) ^ (t) = if" (Vfe_i (t) - 2Vfe (t) + Vfe+i (t)) - Vfe (t) + 3Vfe3 (t) + Ul (t) Vfc (t) + i^2 (a;^, t), fc = l,if-2, 

(5.32) (t) = if" (Vic_2 (t) - 2 Vk-i (t)) - Vic_i (t) + 3V|_i (t) + V|_i (t) Hk-i (t) + i^2 {xk-iU) , 
at 

where (tife (t), Vfc (t)) identically stands for (t), (t)) and the initial conditions are 


(5.33) {eq:5.29}(;7fc (0), Vfc (0)) = (mo (a:fc) ,ho (a;fc)), (Wfc (0), Vfc (0)) = (ui (a;fc) ,01 (a:fe)), fc = 0, if. 

Here we also notice that the values of Uq (t) and V^f (t) are known, so those cases are not considered. 

A basic numerical approach is achieved below by using the linear recursive method where the nonlinear terms are 
linearized. Subsequently, after some rearrangements, we rewrite the linearized differential system to be a differential 
equation where the solution includes all discrete solutions of the linearized system. Such a unifying way also guarantees 
that the numerical solution uniquely exists. So by doing this, at the m-th iterative stage (m > 1) the linearized differential 
system of (I5.27I) - (I5.33|) one by one becomes 


n y('^) 

(5.34) {eq;5.33-l} +3 (t)) V (t)) (t) + Fi (xi, t) , 

it) - 2Ujr^ it) + Uif, (t)) - Ui^'> (t) 

+3 (t))" + (t))" (t) + El {xk,t), k = 2,if-l, 


dU, 


(m) 


dt 


it) = 


(5.35) 


(5.36) 


dU 


(m) 

K 


dt 


{t) = K 


yr"Ht)] -uy 


it) 


- if" ( uy it) - uy it)) -uy it) + 3 it)) + fi ixK,t) , 


.(m) 


























and 


(5.37) 


dV, 


(m) 


dt 


{t) = K 




(m-l) _ y{m) 


+ (vl^'> (t) - (t)) - vlT'’ (t) + 3 {t))' + F2 {xo, t) 


(5.38) 


dt 


it) 


it) - it) + T/i™) (t)) - vf) (t) 

+3 (t))' + (o) (t) + F 2 ixk,t), k = l,K-2, 


,, ,(m) 

-^it) = 7f"(4"!Ui)-24-\w)-vlT-^W 

(5.39) +3 (4^7') (i))' + (u^K-i^ (O)' 44 {t) + F2 ixK-i,t), 

(5.40) {eq:i(0^T^ (0), ^4 (0)) = (^0 ixk) ,io ixk)), (4™’ (0), (0)) = iui ixk),vi ixk )), fc = 0^^. 

Then, a matrix-type differential system derived from (I5.34I) - (I5.40I) can be expressed as follows: 


(5.41) {eq:5.40} (t) = (t) (f) -h it) , it) = K™) (f) (f) -h it) , 

LLL lAj L 

where the solutions and the block matrices x and the functions f4\f4^ € 

are defined by 




T[j(m) 


rr("i) 

K-1 

tA^) 


F^) (f) = 


0 

0 

3 ( 4 '””'^ it)y + Fiixut) 
3 ( 4 ™-') {t)y +F^{x2,t) 


7/(™) 

^K-1 

't/(™) 

*^0 

4™^ 


3 ( 4 ™-/) {t))' + Fi{xK-l,t) 

K (c/4-') (f))' + 3 (c/4-'^ (<))" + F, {xK,t) 

0 

0 


A — 


(m) 


V, 


K-2 

r(m) 

K-1 
■ ) 


V) 


(m) 


Ff) it) = 


77 (4™-') (f)) +3(4’"-'^(f)) +F2(a:o,t) 
3(4’"-'^ it)y + F2ixut) 


,(m) 

K-2 

,(m) 

K-1. 


3(4-4^ it))" +F2ix k- 2, t) 
3(4-7'^ it)) +F 2 ixK-l,t) 























-*) _ 






al := al 




^2 


?,(™) 


(t) = -2K^ + (v, 










(m-l) 


associated with the initial datum 


(m) 

4-1 




-K^ 


'^K-2 


l(™) 

^K-1 


1 0 
0 1 


0 • 

-1 0 

0 -1 


1 0 
0 1 


0 


-1 


0 


0 -1-K 

... 0 


-1-K 0 
0 -1 


-1 0 
0 -1 


^(-) ^ _2^2 ^ 


(m-l) 

k 


k = l,K-l, 


(5.42) ■[eq;5.41} 


( 0 ) = 


uo(a;i) 


f'o (a;o) 

+ (2^2) 


Vo {xi) 

Uo (xK-l) 


Vo {XK-2) 

Uo (xk) 

, ( 0 ) = 

Vo (xk-i) 

Ui (xi) 

Vi (xo) 

Ul (X2) 


Vl (xi) 

Ui (xK-l) 


Vi {xK-2) 

ui (xk) 


yi {xk-i)_ 


Here the initial guess for our linearization method is simply initial values, i.e. (11^°^ (t„), (tn)) = (0), (0)) 

for all n = 0,N. In general, it is remarkable that the system (I5.41I) - II5.42I) might be stiff, then it is not good for our 
implementation since the objective is considering the system in a large time. Understanding the basic instability coming 
from stiff systems, we therefore apply the well-known implicit Euler method which reads 


u: 


(m) 

n+1 


4m) 

'n+1 


= (l- ^A(-) (t„+i)) ' , 

(Wi)) '(^V+) + ^F«(Wi) 



















Figure 5.1. Approximation solutions. 



Figure 5 . 2 . Exact solutions. 


for n = 0, N — 1, where I stands for the 2K-hy-2K identity matrix, and together with the conditions (15.421) . 

Choosing m = 5,T = 20, and K — N = 50, we plot in Fig. 1 the approximation function {u {x,t) ,v {x,t)) solutions 
to our problem (P) considered in this subsection. We surely check that such functions not only behave like the exact 
solutions (j5.24|l (decay exponentially in time), but also completely have the same shapes corresponding to each exact 
solution. Indeed, let us observe those which are sketched in Fig. 2. 

Furthermore, numerical results of the solutions {u,v) together with the exact solutions {uex,Vex) at nodes for 

n = 10; 20; 30, and various values of error in the entry-wise norm 

£n,k{u)= max max \uex {xk,tn) - u{xk,tn)\ , 

l<k<K l<n<N 

= max max \vex ixk,tn) - v {xk,tn)\, 

l<k<K 0 <n<N-l 


£n,k (f) 












n 

ttex ( 1 ) 

u (h^n) 


ttex ( 5 5 ^ ( 5 ’ ^n) 

r 

10 

1.54436330e-03 

2.91855517e-03 

1.37419186e-03 

20 

2.82860006e-05 

7.20712002e-05 

4.37851996e-05 

30 

5.18076174e-07 

1.77972692e-06 

1.26165074e-06 

n 

'^ex 

V (t,in) 


'^ex ( 5 1 tn) ^ { 5 ? ^n') 

1 

10 

3.86090827e-04 

7.29514168e-04 

3.43423340e-04 

20 

7.07150017e-06 

1.80147701e-05 

1.09432699e-05 

30 

1.29519043e-07 

6.22799676e-07 

4.93280633e-07 


Table 1. Numerical results at nodes (|,tn) for n = 10; 20; 30. 


K 

N 

£n,K (u) 

£n,K (v) 

50 

50 

6.68545424e-03 

6.68150701e-03 

100 

100 

3.59475057e-03 

3.59201931e-03 

150 

150 

2.45841870e-03 

2.45632948e-03 

200 

200 

1.86793338e-03 

1.86628504e-03 


Table 2. Numerical results for the loo norm error £n,k- 


are all given in Table 1-Table 2, respectively. To demonstrate the fact that the loo norm error £m,k decreases (and 
obviously tends to zero) as K^N increase without any attention to the discretization (like CFL conditions), we show in 
Table 2 the error values when such constants go from 50 to 200. As expected from analysis, our numerical methods used 
above are reasonable and efficient. 
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